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Abstract 

We have found an effective method of calculating the Wigner function, being a quantum analogue 
of joint probability distribution of position and momentum, for bound states of nonrelativistic 
hydrogen atom. The formal similarity between the eigenfunctions of nonrelativistic hydrogen atom 
in the momentum representation and Klein-Gordon propagators has allowed the calculation of the 
\ Wigner function for an arbitrary bound state of the hydrogen atom. These Wigner functions for 

' some low lying states are depicted and discussed. 
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There are many ways to visualize quantum states of a single particle. Wave functions 
in position representation are a natural tool to present probability distribution in the con- 
figuration space. On the other hand it is easier to visualize momentum distribution if the 
momentum representation of the wave function is used. Phase space description, namely 
Wigner function of a state, provides a natural generalization of joint position and momen- 
tum distribution. Wigner function corresponding to a state described by wave function ip{r) 
is defined as follows [l|: 



(27T)' 



i>*(k + q/2)e~ iqx ^(k-q/2) (1) 



As it is well known the Wigner function for most quantum states is non-positive and bounded 
to an interval -\[— 1, 1], with marginals in f and k corresponding to momentum (wave- vector) 
and position quantum probability distributions. In the classical limit the Wigner function 
becomes a classical phase space distribution. These and other properties of the Wigner 
function with its applications in various branches of physics have been reviewed in a number 
of articles and books 0. 

It is well known that Wigner function given by Eq. can be easily calculated, ana- 
lytically or numerically, for most of one dimensional systems. In case of three dimensional 
problems, especially the ones with spherical symmetry, calculations are usually much more 
difficult. Integrals become quite cumbersome and in most cases impossible to handle analyt- 
ically. Despite the existence of analytical expressions for the hydrogen atom wave functions 
in position and momentum representations the form of the phase space Wigner func- 
tion is unknown. Analytical formula for the Wigner function is not even known for the Is 
state of the hydrogen atom. In the literature one can find only a limited number of papers 
devoted to this subject 0,0, and some of the published results have been achieved using 
some kind of approximate methods that have generated controversies @ . In a different con- 
text, the hydrogen atom has been investigated recently using the Kirkwood-Rihaczek phase 
space representation, which is easy because it involves only products of the momentum and 
position wave functions with a proper phase Q. 

Due to the discovery of quantum phase space tomography, Wigner functions have been 
experimentally reconstructed for quantum states of light, vibrational modes of molecules, 
and superpositions of diffracted cold atoms by a double slit (If. It has been discovered 
recently that the phase space plots of the Wigner function provide a unique visualization 
of the quantum state that can unravel such unique quantum properties like entanglement 
of correlated systems |9f or the phase sub-Planck structures of quantum interference 
Because of all these reasons an analytical formula for the phase space Wigner function, of 
such a fundamental system like the hydrogen atom, can be useful for quantum tomography, 
quantum state diagnostic and phase space visualization of negative structures of quantum 
interference. 

This rather vexing situation regarding the analytical form of Wigner functions for the 
hydrogen atom indicates that the calculation of the phase space representation for the hy- 
drogen requires a new approach based on a new method or a calculational trick to overcome 
the old difficulties. It is the purpose of this Letter to present an analytical computational 
scheme for the calculation of the phase space Wigner function for arbitrary bound-energy 
eigenfunction of the hydrogen atoms in terms of a set of generating differential operators 
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acting on a simple single integral. This general "hydrogen atom integral" (HAI), dependent 
on the Bohr radius and other geometrical parameters of the hydrogen eigenfunction, can be 
easily calculated numerically. 

Before we present the outline of the technical features of the general approach, we first 
summarize our main results, and illustrate the power of our methods showing for the first 
time exact phase space plots of the hydrogen in the Wigner representation. The main result 
of our paper can be written in the form of the following formula for the phase space Wigner 
function for the hydrogen energy eigenvectors: 



W, 
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I(r,k,rk,b 1 ,b 2 ) 



b\=b2=l/na 



(2) 



In this formula D n ; m is a linear differential operator reproducing the Wigner function for all 
hydrogen states from a simple hydrogen atom integral (HAI) defined below by Eq. The 
HAI depends on three scalars only: r = |rj, k = \k\ and fk = rk cos 9. The two arbitrary 
running parameters &i and b 2 , are determined at the end of the calculations only by 1/na, 
where n is the principal quantum number and a is the Bohr radius. The form of this HAI is 



I(r,k,rk,b 1 ,b 2 ) 



i _^ 1 

du exp(4iurk) . 



exp(-2rC(w)) 



where 



C(u) = \jub\ + (1 - u)b\ + 4w(l - u)k 2 



(3) 
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The formula is the central result of our paper. It can be used to generate the Wigner function 
for an arbitrary hydrogen energy eigenfunction. The HAI plays a role of a generating function 
for all Wigner functions of the hydrogen atom. Differential operators D n / m acting on HAI 
give Wigner functions for all bound states of hydrogen atom. 

Before we explain how this central result has been obtained, we write down the formula 
for the Wigner function of the ground state. For the Is state this operator reads: 
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ir 3 a 3 d(bi)d(b 2 2 ) 
and the corresponding Wigner function (J2J is 
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bi=b 2 =l/a 



(6) 



We can perform all the derivatives and write the Wigner function as a single integral, but 
this lengthy formula is rather useless for this letter, and we omit it. We note that the 
Wigner function for the ground state depends only on three scalars r, k and 9. A very 
simple numerical calculations of the HAI leads to the Wigner function for the ground state. 

In Figure (JTJ), we have depicted contours of 4irr 2 k 2 Wioo(r,k,9) for selected values of 9. 
These figures should be compared with the only published numerical results obtained 22 
years ago in reference Q. In the following Figure (j2J) we have depicted the same function 
Wioo(r, k, 9) multiplied by factor r 2 k 2 for 9 = and 9 — |. We see explicitly regions of the 
phase space where the Wigner function is non-positive. 
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FIG. 1: Is state: A7rr 2 k 2 Wi s (f, k). Contour plots to be compared with those from reference 
Dashed lines denote a zero-level, dotted lines denote negative values, separate distance between 
the contours are chosen as in 



9 = 



9 = tt/2 




FIG. 2: The Wigner function of Is state multiplied by factor r 2 k 2 . Only the cross-section for 9 = ^ 
is a positive function, for all others values of the angle between position and momentum vectors 
cross-sections have negative values. For 9 = we see explicitly an oscillating structure. 



The power of our method can be exhibited further if higher states of the hydrogen atom 
are considered. Below we just quote the corresponding differential operators for the 2s state: 
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where after all the calculations we put h = 62 = l/2a. With a little patience or help from 
a symbolic software all the differentiations of the HAI can be performed, and a close form 
expression for the Wigner function W^ 200 (r, k, 6) from Eq. (J2J) can be calculated and plotted. 
In Figure © we present contour plots of the 2s Wigner function for various values of 6. 
Again, the negative values of the Wigner function are clearly seen. 



The next example that we want to present in this paper is the 2p state. In this case 
_ 2exp(2ir£) / 9 1 a n / g x 9 x ( d 2 d \ ^ 
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FIG. 3: The Wigner function of the 2s state. Dashed lines denote a zero-level. Other contours are 
separated by 0.01 

The form of this differential operator indicates that the corresponding Wigner function will 
depend on the scalars (r, k) and two solid angles describing the orientations of r and k. 
The Wigner function of the 2p state is no longer a function of scalars r 2 , p 2 and fk is a 
consequence of the fact that the wave function of this state distinguishes the z axis. As in all 
previous cases a close form expression for the Wigner function W^ 2W (r , k) can be calculated 
and plotted. In Figure (@| we show cross-section of W^ 210 (r, k) for Q\ = 9 2 = and the 
corresponding r 2 k 2 W^ 210 (r, k) plot. 

Finally, in the Figure (0) the Wigner function of 2p state with m = 1 is shown. We 
have plotted the cross-sections of W^ 2U (r, k) for Q\ = 9 2 = | , <fi = and selected values 
of k as a function of r and (p2- These plots show clearly that the maxima of the Wigner 
function are reached for k perpendicular to r, which entirely agrees with classical intuition 
that angular momentum has the maximum value for such geometry. Thus, the semiclassical 
features begin to be visible already in the 2p state. Of course the classical features are much 
more pronounced for larger n. 



In the remaining part of this Letter we illustrate the general calculational scheme showing 
explicitly how one can calculate the Wigner function for the ground state of the hydrogen 
atom. From our derivation it will be clear that the method is general and can be applied 
to all bound states of the hydrogen atoms. Although most of the calculations for higher 
excited states look tedious, simple symbolic differentiation of the fundamental HAI formula 
leads to explicit expression for the the Wigner function with arbitrary quantum numbers. 
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FIG. 4: The Wigner function of 2po state: Plot a) presents the Wigner function for 9\ = 62 
b) shows the same cross-section with the Wigner function is multiplied by r 2 k 2 . 
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FIG. 5: The Wigner function of 2p\ state (we have chosen 9\ = 62 = §, <pi = 0). Red line denotes a 
zero-level. Plots a) and b) show cross-section of the Wigner function for fixed values of momentum 
k = 0.1 and k = 0.2, respectively. On the axes are r and the angle between r and k. It is seen 
that maximum is obtained when position and momentum vectors are perpendicular, which agrees 
with classical intuition that angular momentum has the maximum value for such geometry. Plot c) 
presents similar cross-section of the Wigner function multiplied by r 2 k 2 for k = 0.5. As we would 
expect maximum is obtained for ip = 5 and r ~ 4. 



We have found useful for our calculations to work with wave functions in momentum 
representation |3j]. For the wave function of the ground state we have: 



^100 (fc) 



W 



1 + k 2 a 2 ) 2 



(9) 



It is a regular function that for large k decreases as (ka) 4 . Inserting this expression into 
the definition of the Wigner function in the momentum representation, Eq. (0), we obtain: 



2 9 na 



d 6 q — 



exp(— 2if(q — k)) 



(10) 



The power of q in the denominator can be reduced with the help of differentiation over 
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parameters. The following expression for the Wigner function is obtained: 

Wl (fk)- 2 ° 2 fd 3 a ^t?MdM fin 
* 0o( ' } " d(b 2 )d(b 2 ) J dq (bl + + (q- 2k) 2 ) ' (U) 

where b\ and b<i are running parameters to be fixed by the Bohr radius at the end of all 
calculations. Now comes the key element of the calculation. We recognize, that apart of the 
phase factor, the structure of the integrand in Eq. (jllj) has a remarkable formal similarity to 
the product of two Klein-Gordon propagators of quantum field theory in momentum space 
[Til ]. Due to this analogy, we shall proceed with our calculations using the standard propa- 
gator disentanglement techniques introduced by Feynman and represented by the following 
identity 

If 1 1 

' du 2- ( 12 ) 



AB Jo [uA+(l-u)BY 

with A = b\ + q 2 and B = b\ + (q — 2k) 2 to rearrange the integral in Eq. ([lljh The integral 
becomes: 

2 d 2 f 1 
W fl00 (?,k) = ^3 d (h\)d{bl) Jo dU 



x / d 3 exp[-2i(q-k)r\ 

[u{bl + q 2 ) + (1 - u){bl + (q - 2k) 2 )} 2 ' 

After rearrangement of terms and substitution: s = q — 2(1 — u)k we get 

2e 2i ^ d 2 f 1 
W^ 100 (r,k) = -g-j d{h 2 )d{pl) J o du 

3 exp[-2i(s + 2(1 - u)k)r\ 



J [s 2 + m/3 + (1 - u) 7 + 4(1 - «)ufcf ' 1 J 

We recognize in this expression the function C(u) introduced in Eq. (|4*jl. This simplifies the 
notation and the integral (|14j) then becomes: 

-> 2c 2 ^ r 2 -> f exof 2isv) 

k) = ^rg^^y^ /„ <W-4i(l - y d». (15) 

Fortunately, the integral over d 3 s is elementary, we can use the following formula: 

3 exp(-2zsr) tt 2 

rf s WTcW? = W) exp( " 2rC(M)) - (16) 

As a result of all these steps the only remaining integral is over u. The final formula for the 
Wigner function is thus given by: 

^-(r,^) = -^^p^yX d«exp(<4u**0^exp(-2rC(t*)) 



d d d • 

W^ 100 (f,k) = B 100 [-^,—, — ) I(r,k,fk,b 1 ,b 2 ) 



(17) 

61=62=1/0 



as it has been advertised in Eq. (J2J. 

This method works for arbitrary state of the hydrogen atom. The key concept in such 
calculations is to express a given state with quantum numbers (nlm) in momentum rep- 
resentation as a differential operator acting on the ground state followed by a change of 
scale. It is worth noting that these differential operators form an elegant group theoretical 
structure, explained in e.g. 

In conclusion, we have found and exploit formal similarity between the eigenf unctions 
of nonrelativistic hydrogen atom in the momentum representation and Klein-Gordon 
propagators. This allowed us to find and discussed the Wigner function for arbitrary bound 
state of hydrogen atom. 
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